Fractional flux periodicity of a twisted planar square lattice 
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We present fractional flux periodicity in the ground state of planar systems made of a square 
lattice whose boundary is compacted into a torus. The ground-state energy shows a fractional 
period of the fundamental unit of magnetic flux depending on the twist around the torus axis. 



The Aharonov-Bohm effect pj shows that a single elec- 
tron wave function has a fundamental unit of magnetic 
flux <I>o = 27r/e, where — e is the electron charge. We will 
use the units of % = c = 1. The electric and magnetic 
properties of materials are governed by many electrons, 
and each constituent has the above-mentioned periodic- 
ity. However, the fundamental flux period of a material 
need not be <I>o- For example, superconducting materials 
exhibit a period one half the single-electron flux quan- 
tum, that is, 4>o/2, which can be understood by charge 
doubling due to the Cooper pair formation. This is clear 
if one imagine that a fundamental particle (or quasipar- 
ticle) has a — 2e charge due to the attractive interaction 
and that the fundamental flux period of a material is 
equal to that of the quasiparticle (27r/2e) in the system. 
So, one may ask: Is it possible that fractional flux pe- 
riodicity ($0/2, $o/3, • ■ • ) realizes in the absence of the 
interaction? We answer this question for torus geome- 
tries whose surface is a planar system made of a square 
lattice. To analyze the flux periodicity of the systems, we 
will consider the persistent currents which are suit- 

able to examine the flux periodicity of the ground state 
since they are regarded as the Aharonov-Bohm effect in 
solid state systems 0. 

The quantum mechanical behavior of conducting elec- 
trons on a square lattice is modeled by the nearest- 
neighbor tight-binding Hamiltonian: 
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where t is the hopping integral, A c * is a constant exter- 
nal gauge field (vector potential), and T a (a = x, y) is the 
vector connecting each site in the direction of x and y. 
a,i and a\ are canonical annihilation-creation operators of 
the electrons at site i that satisfy the anti-commutation 
relation {ai,aj} = 5y. The energy eigenvalue of the 
Hamiltonian is parametrized by the Bloch wave vector 
k as 
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The geometry of a torus can be specified by two vec- 
tors: chiral Ch — NT X + MT y and translational T w — 
PT X + QT y vectors (we borrow this terminology from the 
carbon nanotube context 0), where N,M,P, and Q are 
integers. It is useful to rewrite T a in terms of the chiral 
and translational vectors as 
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where we define N s = NQ — MP. In this study we will 
take x axis in the direction of Ch (M — 0), and denote 
P = 5N to describe the twist around the tube axis. Then, 
C h -T w = SNNa 2 where o(= \T X \ = \T y \) is the lattice 
spacing. Figure^Ja) illustrates a twisted torus. It should 
be noted that a torus can be unrolled to a parallelogram 
sheet as depicted in Fig. IHb). In the figure, the two 
lines extending upward from 'u' and downward from 'd' 
and having the same 'x' at the junction are not joined 
for a twisted torus, shown in the lower inset of Fig.^a). 
There are SN square lattices between the two lines. 
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FIG. 1: Twisted torus (a), its net diagram (b) with an exter- 
nal gauge field, and an example of a twisted torus (c). It is 
convenient to consider a parallelogram as the net of a twisted 
torus. We draw several square lattices in (b). 

The Bloch wave vector satisfies the periodic boundary 
condition through which the geometrical (or topological) 
information, such as the twist, is put into the energy 
eigenvalue. We decompose wave vector k as (iiki + iiiki, 
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where fj,± and fj, 2 are integers. Here, k± and k 2 are defined 
by 

C h -k x = 27r, C fc • fc 2 = 0, T w -h = 0, T w • k 2 = 2tt. (4) 



By means of Eqs.© (with M 
we rewrite Eq.© as 



and P = SN) and 0. 



where P2( = ^ik-i) denotes the momentum along the axis. 
Notice that A cx • T TO can be included by replacing p 2 with 
the covariant momentum and Eq.JHJ) is for p 2 > 0. As 
for p 2 < 0, the energy dispersion relation is defined by 
—'H f j, 1 up to 0(p 2 ). The coefficient of p 2 defines the Fermi 
velocity and that of p\ the effective-mass if vf(pi) = 0. 
They are defined respectively as 
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where we assume A cx ■ Ch — 0, and A cx ■ T w corresponds 
to the Aharonov-Bohm flux penetrating the center of the 
torus. It is important to note that the twist works as 
an extra gauge field, ^4 twlst ; along the axis because 
the second term of the right hand side of Eq. (JSJ can be 
rewritten as exp ^i 27TfJ/2 ~( eA — hZk^ i where we de- 
fine A twiBt -C h = and A twist ■ T w = 2ir5N/N. Like A GX , 
^twist can gj^jf^ the- wave vector. However, the shift de- 
pends on the twist (SN/N) and the wave vector around 
the tube axis (fi±ki). Each energy eigenvalue is deter- 
mined by pL\ and p 2l and the conducting electrons with 
an integer p\ (—N/2 < \l\ < N/2 — 1) form an energy 
band. This means that coupling between A twlst and the 
conducting electron preserves the time-reversal symme- 
try of the whole system. This is contrasted with the 
broken time-reversal symmetry by A cx . 

Persistent currents are defined by differentiating the 
ground-state energy (Eq(Q)) with respect to the mag- 
netic flux $ = A ex ■ T w : 
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Eo(3>) consists of the energy eigenvalue of the valence 
electrons and has the period of $o- Therefore, I pc (®) 
also has the period of <!>o and is known to be a saw- 
tooth shape as a function of <!>. However, to calculate the 
persistent currents in long systems \T W \ 3> \Ch\ (or Q 3> 
N), we do not need to sum over the energy eigenvalue 
of all the valence electrons. We can just calculate the 
Fermi velocity of each energy band. In terms of Fermi 
velocity, the amplitude of the persistent current for the 
/ii-th energy band is well approximated by Q 
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where \T W \ is the system length and vf(^i) denotes the 
Fermi velocity of the pi-th energy band. We have ne- 
glected a higher order correction to the amplitude of 
order of I(jii)0(N/Q). We fix fi\ and expand Eq.© 
around the Fermi level (we assume half-filling and set the 
Fermi level as E-p — 0) to obtain the energy dispersion 
relation of the /iith energy band as 
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We now consider the persistent currents in an un- 
twisted torus (SN — 0). In this case, all energy bands 
yield the same function for the persistent current, that 
is, a saw-tooth curve as a function of <E> with the same 
zero points (the flux for which the amplitude of the cur- 
rent vanishes) but with different amplitudes. In fact, the 
amplitude of the total current is given by a summation 
of all amplitudes: 
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Persistent current / pc in the torus is given by Tot mul- 
tiplied by <f>/ir (cf> = 27r(<f>/$ ))- The linear relation be- 
tween $ and the persistent currents has a periodicity of 
$0; and we then have a saw-tooth current: 
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The saw-tooth persistent current is a characteristic of 
the non-interacting theories at zero temperature. Each 
saw-tooth curve loses its sharpness due to disorder, or at 
finite temperature Q. 

Next, we analyze a twisted torus (SN ^ 0, see 
Fig. HJc)). As we have already shown, the twist behaves 
as an extra gauge field and shifts the wave vector along 
the axis direction, so that the zero points of the persis- 
tent current also shift 0. As a result, in order to calcu- 
late the total current, we must sum the saw-tooth curves 
that have different zero points and different amplitudes, 
both of which depend on p,\. This can be achieved in a 
straightforward manner, and the formula for the persis- 
tent currents in a twisted torus is given by 
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(8) where Iff ((f)) is expressed as 
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Persistent currents exhibit the following characteristics 
depending on an even or odd number of SN. Term 1 + 
cos(iTnSN) in Ea. l|14|) vanishes if nSN is an odd number 
or is equal to 2 for other cases. Then Ea. (|14f> reduces to 



Finally, let us remark that the amplitude of the current 
exhibits a nontrivial dependence on $ when N — > oo with 
a fixed value of SN. If we divide Ea. l|15|l by N and take 
the limit of N — > oo, we then have 
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where we have introduced the sequences 
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The difference in periodicity for cj> is clear because 
sin(2n</>) has a period one half the fundamental unit of 
magnetic flux, or $o/2- This conclusion is valid for any 
number of N, but a fundamental flux periodicity for a 
current demonstrates more interesting behavior for a spe- 
cific combination of N and SN. 

We consider a particular structure for 8N/N =1/3 
where SN is an even number. In this case, the funda- 
mental period of C e n is 3, i.e., C* +3 = C*. When N > 1, 
we have C{ = Cf « -8tt/37V and Cf « 2^/tt. We then 
obtain 
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where in the last line of Ea. l|18J) we assume — 7r/3 < <fi < 
7r/3. The fundamental flux period of Ea. (|18|l becomes 
$o/3- It should be noted that this argument can be ap- 
plied to other twisted structures. For instance, when 
SN/N = l/Z (l/Z = 1/4, 1/5, •■•), we obtain another 
fractional period of $o/Z. To our knowledge, this is the 
first time a fractional periodicity for a twisted boundary 
condition has been found. We note that the total persis- 
tent currents are still saw-tooth curves as expected from 
the non-interacting theories. 
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When SN ^> 1, we sum n in the above equations and 
obtain 
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where we have used the following mathematical formula: 
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In Fig. [3 we plot the persistent currents of Eq.J^OJ) for an 
even and odd number of SN as a function of </>. It should 
be noted that the persistent currents are not standard 
saw-tooth curves even though we are considering non- 
interacting electrons, and the functional shape of Eq.J20|) 
for SN is an odd number similar to sin(2</>) but not iden- 
tical (see Fig. EJ) . 

Let us summarize our results: (1) A fundamental flux 
period of the ground state can be generally fractional as 
§o/Z, depending on the ratio of twist SN to N, although 
we do not assume any interaction that can form a quasi- 
particle of charge Ze; (2) As N — > oo with fixed value 
SN(~^> 1), the currents are not standard saw-tooth as is 
normally expected from non-interacting theories; and (3) 
To observe fractional periodicity it is essential that the 
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FIG. 2: Persistent currents as N — > oo f Eg 1201 for an even 
(solid line) and an odd (dashed line) number of 5N. We 
multiply 2 (16) for an even (odd) number case. We also plot 
0.8 sin(20) (dotted line) for comparison. 



system contain many electrons: N 3> 1 and Q ^> N, 
indicating a kind of "many body" effect. 

The primary factor in these results is Fermi surface 
structure of the square lattice. It appears that many 
energy bands cross the Fermi level when we roll a planar 
sheet into a cylinder, and electrons near the Fermi level 
in each band contribute to the persistent currents, which 
interfere with one another due to the twist. This point 
can be made clear by comparing a square lattice with 
a honeycomb lattice (which possesses only two distinct 
Fermi points). In this case, we can observe at most a 
4> /2 periodicity @|. 

Here, we discuss some possible extensions of our re- 
sults. Because we have observed a fractional periodic- 
ity in the ground state, one may ask the following ques- 
tion: "Is it possible that the persistent currents show 
multiple periods such as 2<E>o or 3<&o as a fundamental 
period?" To answer this question, we consider higher 
genus materials (g; number of holes), the ground state 
of which exhibito<I>o periodicity depending on the genus 
(g = 2,3, • • • ) 8j. A planar system comprising a finite 
number of square lattices is an example of higher genus 



material. We note that the conducting electrons in that 
case are also assumed to be non-interacting and have 
a single electron charge. Moreover, other examples are 
known where the charge of a quasiparticle itself becomes 
fractional. The quasiparticle in the quantum Hall effect 
has fractional charge e/3 |j and there is a model con- 
taining a fractional e/2 charged soliton in 1+1 dimen- 
sions [l0|. Those systems might exhibit 3$o or 2<E>o pe- 
riodicity in the ground state respectively. However, frac- 
tional periods other than $o/2 (which may correspond 
to a multiple charge e — + 3e, 4e, • • • ) have thus far not 
been observed as far as we know. 

So far we have ignored the dependence of persistent 
currents on the value of Q and implicitly assumed that 
Q is a multiple of N. For a general Q, the persistent 
currents are not so simple |6j as the results obtained in 
this study. For example, when the reminder oiQ/N is an 
odd number, the one-half periodicity may emerge even in 
the absence of twist, moreover, another fractional period 
may appear for a specific value of Q. 

In summary, we have examined the flux periodicity of 
the ground state of the conducting electron on a torus 
of square lattice. We found that the persistent currents 
and the ground-state energy of the systems show frac- 
tional periods of the fundamental unit of magnetic flux 
($o/2, $o/3, • • • ) depending on twist SN and N. Fur- 
thermore, for the case of N — * oo with a fixed value of 
SN(^> 1), the persistent currents are not standard saw- 
tooth curves as expected from the non-interacting theo- 
ries at zero temperature. 
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